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The standard system-plus-reservoir approach used in the study of dissipative systems can be mean- 
ingfully generalized to a dissipative coupling involving the momentum, instead of the coordinate: 
the corresponding equation of motion differs from the Langevin equation, so this is called anoma- 
lous dissipation. It occurs for systems where such coupling can indeed be derived from the physical 
analysis of the degrees of freedom which can be treated as a dissipation bath. Starting from the 
influence functional corresponding to anomalous dissipation, it is shown how to derive the effective 
classical potential that gives the quantum thermal averages for the dissipative system in terms of 
classical-like calculations; the generalization to many degrees of freedom is given. The formalism 
is applied to a single particle in a double-well and to the discrete </> 4 model. At variance with the 
standard case, the fluctuations of the coordinate are enhanced by anomalous dissipative coupling. 



I. INTRODUCTION 

The system-plus-reservoir (SPR) model is the 

most common and effective approach to the treatment 
of quantum dissipation. This model realistically assumes 
the dissipation in a system as due to the interaction 
with a reservoir (or bath, or environment): anyone of 
the numerous degrees of freedom of the reservoir is only 
weakly perturbed, so the reservoir is at thermal equilib- 
rium; moreover, in order to reproduce the general phe- 
nomenological dynamics (quantum Langevin equation for 
the coordinate) it is sufficient to assume a linear inter- 
action and a bath of harmonic oscillators [Q , so that in 
the imaginary-time path-integral formalism it is possi- 
ble to integrate out the bath variables and get a reduced 
description of the system in terms of a bilocal influence 
action . Starting from this framework, in previous pa- 
pers we obtained the classical effective potential for the 
calculation of thermal averages in a dissipative quantum 
system Q and gave applications to the 4 -chain model 
and to a 2D array of Josephson junctions 

At variance with the standard approach of coupling 
the bath with the system's coordinate, we consider in 
this paper the complementary possibility of a coupling 
with the momentum, that Leggett called anomalous dis- 
sipation §. In this case the dynamics can be reduced to 
a pseudo-Langevin equation where the dissipative term 
contains the second derivative of the potential. This 
raises the question whether one can tell the dissipation 
mechanism just from the phenomenological behavior: we 
believe that one should have a physically meaningful mi- 
croscopic model for the bath rather than rely upon the 
phenomenological counterpart. The case of anomalous 
dissipation occurs, for instance, when the effect of the 
blackbody electromagnetic field on a Josephson junction 
is considered fiof] . 

As for thermodynamics, we will show that the anoma- 
lous influence action can be derived following the original 



Feynman- Vernon pd[ idea, which, due to the momentum- 
path dependence, can be pursued at the price of involving 
the full phase-space path integral. 

The main goal of this paper is in the treatment 
of this path integral within the effective-potential for- 
malism [Q, that reduces the evaluation of quantum- 
dissipative thermal averages to much simpler classical- 
like configuration integrals. This is accomplished also 
in the case of many degrees of freedom, making quanti- 
tative calculations feasible. In general, the inclusion of 
dissipation by coupling the system's momenta with the 
environment results in higher quantum fluctuations of the 
coordinates, while those of the momenta are suppressed. 

We introduce in Section [n] the concept of anomalous 
dissipation; it turns out that the influence action depends 
on the momentum path and it is hence nec essa ry to use 
the Hamiltonian path integral. In Section III we treat 



it within the pure-quantum self-consistent harmonic ap- 
proximation |12[ , deriving the classical effective poten- 
tial and the corresponding classical-like expressions for 
quantum thermal averages; in the case of many degrees 
of freedom the treatment is simplified by the additional 
'low-coupling' approximation. Eventually, in section IV 
we show how the framework works for the single particle 
in a double-well potential and for the </> 4 -chain. 



II. FROM STANDARD TO ANOMALOUS 
DISSIPATION 



The standard SPR Hamiltonian has the form 

-mi 
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where p and q are the momentum and coordinate of the 
'system' (here, one single particle in the one-dimensional 
potential V), while pi and qi are those of the ^-th de- 
gree of freedom of the reservoir. It is quite a general 



result Hi that a suitable distribution of the bath param- 
eters {m,£,uj£} (see Appendix [A]) can reproduce, when 
the bath variables are eliminated from the equations of 
motion, the most general quantum Langevin equation 
namely, 

m'i + f dt' j(t-t') 'q(t') + V'(q) = F(t) . (2) 



Since it does not explicitly contain the 'microscopic' de- 
grees of freedom, the Langevin equation constitutes a 
macroscopic description of the dissipative system with a 
clear phcnomcnological meaning, in the sense that the 
memory function ^y(t) can be thought to be determined 
experimentally. The thermodynamic density matrix for 
the standard SPR model at the equilibrium temperature 
T — has the path-integral expression (A.2), where 
dissipation is described by the additional influence ac- 
tion 

f m du f ph du' 

sA<i( u )} = y o 2K J ~Jm k ( u ~ u } i( u )i( u ') ' ( 3 ) 

and the Matsubara transform of the kernel k(u) is di- 
rectly related with the Laplace transform of the memory 
function as k n = \v n \ j(\v n \), where v n = 2im/ (ih. For a 
harmonic potential one obtains for the density matrix the 
exact result (A5), so it appears that standard dissipation 



quenches (q 2 ) and rises (p 2 ). This breaks the canonical 
symmetry between the coordinate and the momentum, 
and is obviously a consequence of the bath being coupled 
to the coordinate in the standard SPR model (|l|). In the 
general case of a nonlinear int eraction V{q) the evalu- 
ation of the path integral (|A2|) for the SPR model be- 
comes quite complicated, and was the subject of Refs. || 
and 0, where a classical effective potential suitable to 
reduce the problem to classical-like expressions was in- 
troduced as a generalization of the original approach of 
Refs. pi and p). 



The different case of anomalous dissipation occurs 
when the SPR model is modified including the momen- 
tum in the place of the coordinate in the interaction with 
the reservoir. The dependence upon the bath coordinates 
and momenta is quadratic, so they can be exchanged at 
will and the anomalous SPR Hamiltonian can be written 



H = ^ + V(q) 
Im 



(4) 



It is immediately seen that Eq. (22) of Ref . |]To[ , regarding 
a single Josephson junction interacting with the black- 
body electromagnetic field (in the dipole approximation), 
can be cast exactly in this form, constituting a first phys- 
ical example of the relevance of treating anomalous dis- 
sipation. 

Let us look at the dynamics of the system (||): as in 
the standard SPR case, one can derive the equations of 
motion and then eliminate the bath variables. The result 



is a pseudo-Langevin equation, 



+ m dt' V (t-t') dtiV'(q{t')) + V'(q) = F(t) . 
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To interpret Eq. (||) , observe that in the classical limit one 
has dfV'{q{t')) = V" (<?(*')) #(*')> so that in thc neigh- 
borhood of the stable minimum, where V"(q) ~ mu: 2 
tends to a positive constant, it reduces to the standard 
form of the Langevin equation (|^); for a harmonic po- 
tential this always occurs. The physical difference is that 



the full damping function ^{t) 



2 i](t) depends on 



the system's mass and potential, not barely on the reser- 
voir's characteristics, so that the frictional force depends 
on how the system's motion arises, which prevents from a 
simple phenomenological interpretation JlJ . On the same 
footing, the concept of anomalous dissipation was consid- 
ered by Caldeira and Leggett (see Appendix C of Ref. H ) 
and further analyzed by Leggett |9[. Their purpose was 
to establish that in the context of tunneling problems 
the standard SPR model (Q) (with a possible coupling to 
a nonlinear function of the coordinate) is the most gen- 
eral one to be considered in order to infer the effects of 
damping from a knowledge of the phenomenological qua- 
siclassical dissipative equation (|); therefore, they rule 
out the anomalous case as it can give a negative friction 
coefficient (V" < 0) in the pseudo-Langevin equation 
over much of the tunneling region, which would lead to 
qualitatively different results. On the other hand, since 
the pseudo-Langevin equation reduces to the Langevin 
one in the dynamical asymptotic region (i.e., close to 
the equilibrium configuration), a purely phenomenologi- 
cal approach to dissipation based on the dynamical linear 
response cannot distinguish whether the underlying dissi- 
pative mechanism is of the standard or of the anomalous 
type, while this distinction plays a fundamental role in 
the quantum statistical mechanics. However, rather than 
being a schematization for a phenomenological dissipa- 
tive behavior, the mechanism can happen to be deter- 
mined from a physically sound microscopic model which 
can have the form ([!]) : once such microscopic SPR Hamil- 
tonian is known, there is no point in trying to reduce 
it to a phenomenological description. Although exam- 
ples can be given of how anomalous dissipation can arise 
just reinterpreting the canonical variables in a standard 
SPR system (e.g., a linear dissipative system as a simple 
RLC circuit P|, using the canonical symmetry between 
magnetic flux and charge), it is more significant that the 
anomalous dissipative mechanism is obtained on a mi- 
croscopic basis for Josephson junctions, when quantum 
electromagnetic fluctuations are considered in the dipole 
approximation Jl0| . 

In order to get a first insight about the dissipation ef- 
fect onto the thermodynamics, let us consider the har- 
monic oscillator. In Appendix |b| the analytic solu- 
tion is easily obtained using the canonical symmetry 
between coordinate and momentum. At variance with 
the standard case, it appears that anomalous dissip atio n 
quenches the fluctuation of the momentum, Eq. dB^), 
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and enhances those of the coordinate, Eq (|B10|), as it 
is expected since the environment now performs a 'mea- 
surement' of the momentum. This kind of behavior is 
expected also when the system is subjected to a non- 
linear potential ||, but in this case the evaluation of 
quantum-dissipative thermal averages is by no means 
straightforward, and one has to resort to some approxi- 
mation scheme. The purpose of the next Section is indeed 
to derive the effective classical potential and the related 
prescriptions in order to evaluate the system's thermo- 
dynamic averages for the anomalous SPR model (f|). 

In writing the path integral for the density matrix cor- 
responding to the Hamiltonian (^), it appears that in 
order to integrate out the bath variables one must re- 
sort to the full phase-space formulation for the system, 
so that the momentum path p(u) enters the interaction 
with the bath. Indeed, the anomalous influence action 
turns out to be (see Appendix ^) 

SMu)} = k{u-u>) p(u)p(u>) , (6) 



o 2n. J 



and the kernel k(u) = K n e tL ' nU is now directly related 
with the Laplace transform of the memory function rj(t) 
as 



v(Wn\) 



(7) 



Note that k(u) is periodic and, since K n — n- n , also even, 
so one has the properties k(u) = k{—u) = k{u + 0h); in 
addition kq = 0, meaning that the influence action does 
not include a local component, which trivially renormal- 
izes the mass. 

For the phase-space path integral it is convenient to 
use the symmetric (Weyl) ordering prescription [|l5|Jl6| 1. 
The Weyl symbol for the density operator can therefore 
be written as 



p(p,q)= / V\p u ,q a ] e-Sb^l-W 



(8) 



where for compactness the arguments are indicated as 
subscripts, and the nondissipative part of the system's 
action reads 



S[pu,q u ] 



0h du 

h 



(quPu-Puqu) 



Pu 

2m 



V(q u 



+-^(qoPp-PQqp) + -^[(q0-qo)p-(p[3-po)q\ ■ (9) 

Note that the arguments p and q do appear explicitly in 
the action, while there are no constraints upon the paths; 
in particular, the end points (po, pp, qo, qp) are also in- 
tegrated over HU . The thermal average of an observable 
O can be written in terms of its Weyl symbol, 



0(p,q)= f dx e -*P*/^ g +||o| g -|) , 



a.s 



(10) 



(11) 



III. EFFECTIVE POTENTIAL 

A. One degree of freedom 

The strategy [[l2| to approximate the density (g) is to 
use the quadratic trial action So[p u , q u ; q] expressed as 
Eq. ([}]) where V(q u ) is replaced by the trial 'potential' 

Vb(g«5 q) = w{q) + - muj 2 (q) (q u - qf , (12) 

which is indeed a functional through its dependence on 
the path's average point 



q[q u 



o 



(13) 



and to optimize its parameters w(q) and uj 2 (q). 
The reduced density 

po(p, q; q) = J v\pu, qu] s(q- q u ) 

x e -So]pu,quiQ]Sj\}hi] fi£\ 

collects the contribution of classes of paths that share 
the same average point q, so the path integral (|^) for the 
trial action So can be cast into the form 



Pa(p,q) = dq p (p,q;q) ■ 



(15) 



What we are going to do is to optimize the reduced den- 
sity p. This is an important point, that justifies the ac- 
curacy of the approximation: we choose the best approx- 
imating Gaussian distribution for the separate contribu- 
tion of each class of paths (labeled by q) , which describes 
the purely quantum-dissipative contribution to the fluc- 
tuations around any 'classical' configuration q. In other 
words, the classical fluctuations are accounted for ex- 
actly. 

The reduced density can be evaluated analytically: by 
including the constraint for the variable q in the action 
through the Fourier representation of the Dirac delta 
function it is possible to get an integral expression in 
terms of the dissipative harmonic oscillator density ma- 
trix. In Appendix |fj we report the actual calculations 
that lead to the final result, which is a Gaussian distri- 
bution in (p, q): 



po(p, q; q) 



V/2A p -(q-q) 2 /2a 



where the effective potential reads 



1 



V eS (q) = w(q) + -qV{i) i 



with 



n=l 



v\ + (1 + mn n )uj 2 (q) 



(16) 



(17) 



(18) 
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and the relevant variances are 



\(a) = - V ^'(g) 



1 + WiKi 



^ ^ ^ + (! + mK n) u 2 (q) ' 



(19) 
(20) 



apart from the dependence on the variational parame- 
ter Lu 2 (q), thes e var iances have the harmonic-oscillator 
form (BE) and ( B10| ) with the important difference that 
a lacks the n—0 term, i.e., its classical contribution. 

We introduce henceforth the coordinate fluctuation 
variable £ = q — q in the place of q and the following 
double-bracket notation for the Gaussian averages over 
the variables p and £ defined by po : 



(p 2 )) = X(q) , (O) = a(q) 



(21) 



The physical interpretation of the above formulas be- 
comes transparent if one expresses the resulting approx- 
imation for the thermal average ([ll]) using Eqs. ( fijl ) 
and (|l6|), namely 



understood as not operating onto its coefficient, is such 
that AF(q) = (F(g+£))), and allows us to rewrite the 
effective potential in a more compact way: 

V og (q) = (1 - A) e A V(q) + - . (27) 

P 

It appears that the first term gives corrections of order a 2 
to the potential, while first order renormalizations arise 
from the last term. 



B. Many degrees of freedom 

Let us now consider a general system with N degrees 
of freedom, i.e., canonical coordinate and momentum op- 
erators q — {cji}i=i n and p = {pi}j=i,...,jv, with the 

commutation relations [<?i,Pj] = iSij (we set ft = 1 in 
this subsection), and described by a Hamiltonian with a 
quadratic kinetic energy and a nonlinear potential term, 



H = \ t pA 2 p 



V(q) ■ 



(28) 



{0) ~ Z V 2nh 2 (3 



dq e 



(o(p,q + 0) ; 



(22) 



the average over £ accounts for the nonclassical fluctu- 
ations of the coordinate, while the classical part is ac- 
counted for exactly by the classical-like expression for the 
thermodynamic average; the average overp keeps into ac- 
count the full fluctuation, since its classical contribution 
is exactly Gaussian. The detailed discussion made in 
Rcf. applies in this case as well, provided one keeps 
into account that the contribution of the dissipative non- 
local action is also accounted for by the variances A and 
a. 

The last point concerns the optimization of the param- 
eters w(q) and Lo 2 (q). This can be performed imposing 
that the trial and the true potential (and their derivatives 
up to the second one) have the same p -average 



i 



(V(q+0) = w{q) + -mij 2 (q) a(q) , 
(V"(q+0))=mcj 2 (q) , 



(23) 
(24) 



while the condition for the first derivative is trivial, since 
a linear term ~ (q—q) in Vq does not contribute to the ac- 
tion by the very definition of the average point. Eq. ( p4| ) 
is a self-consistent equation that identifies the parame- 
ter u) 2 (q), while the determination of w(q) in Eq. ( f23|) 
permits to rewrite the effective potential as 



The differential operator 



A=±a(q)d, 



(25) 



(26) 



The real matrix A = {^4fj} is symmetric, *A = A, and 
positive definite. The influence action for anomalous dis- 
sipation takes the form 



f^dn f 13 dn' 

sAp(u)} = J o yJ T tp{u) K{u - U ' ] p{u ' ] ' 



(29) 



where the kernel matrix K(u) = {^(w)} is a real sym- 
metric matrix that replaces the scalar kernel n(u) of the 
single-particle function of u it keeps its sym- 

metry and periodicity, K{u) = K{—u) = K(f3—u), and 
satisfies duK(u) = 0. For instance, in the case of N 
independent identical baths coupled to each momentum 
Pi one simply has a diagonal kernel, Kij(u) — Sij n(u); 
this cas e will be considered for the application shown in 
Section IfVC. 



The Weyl symbol for system's density matrix is ex- 
pressed as an TV-dimensional integral, 



P{v,<l) = J dq p{p,q;q) 



(30) 



in terms of the reduced density collecting the contribu- 
tions of paths sharing the average configuration q, 

x e - 5 [p u ><? u ;p.<3]- s [p ll < < J u ;9]- s i[Pu] ■ (3i) 

here the action is separated in a part containing the ex- 
ternal variables (p, q) , 



S[p u ,q u ;p,q] = 



du 



2 -K<Z«) 



AoP/j-VoC/s) + * [Xqp-Qo)P - XPp-Po) q] > (32) 
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and the part containing the system's Hamiltonian, 



S[p u ,q u ;q} = 



du 



^ t PuA 2 p u + V(q u ) 



(33) 



In order to evaluate the effective potential, the trial 
action Sq is defined as S by replacing V(q u ) by 



V (q u ; q) = w(q) + - \q u -q)B 2 {q){q u -q) 



(34) 



with a real symmetric matrix B 2 (q). The calculation of 
the corresponding po , reported in Appendix ||, is not a 
trivial extension of that for one degree of freedom, but 
can take advantage of the results obtained in Ref. for 
the standard SPR model. 

The simplest way for writing the final result is to give 
the expression of the thermal average of a generic ob- 
servable O in terms of its Weyl symbol 0(p,q). This 
fundamental formula approximates quantum averages by 
means of a classical-like expression with the effective po- 
tential V c ff, 



1 



Z\2nf3j detA 



dq e^ v ^(O(p,q+0) , 
(35) 



where (( • )) is the Gaussian average over the varia bles p 
and £ determined by po, as reported in Eq. ( E10 ), and 
can be uniquely defined through its moments 



(et} = c(-q) = W B - 1 *" 



(36) 



(37) 



whose components are the renormalization coefficients 
and 



* n (q) = B (A 2 + K n ) B 
The effective potential reads 



V eS (q) = w(q) + - p(q) 



ith 



^ det(ijj + gn) 
p{q) = ^ In 



n=l 



(38) 



(39) 



(40) 



To determine the parameters w and B(q) we require, as 
in Eqs. ( |23| ) and (pi|), that the parameters of the trial 
action are such to match the po-averages of the original 
and the trial potential, and the same for their second 
derivatives: 

{V{q + i)))=w{q) + \^[B 2 {q)C{q)] , (41) 



((d q AiV(q + S)}=B* j (q) 



(42) 



The latter equation together with Eq. ( Pq) self-consis- 
tently determines the solution for the matrices B(q) and 
C(q). Its matrix character makes it useful to intro- 
duce the 'low-coupling' approximation (LCA), in the very 
same way of Ref. JjJ , so it is sufficient to write here the 
final results. If q Q is the configuration that minimizes 
V e s(q), one has the LCA effective potential 



V eS (q) = e A V(q) - Ae A V(q a ) + /T 1 p 
where p = p(qo) and the operator 



(43) 



(44) 



with C = C(q ), is such that e A V{q) = (Vj(o+£)) 
(within the LCA). Therefore Eqs. (|3q ) and ( |42| ) have 
to be solved only for the minimum configuration, with 
a great simplification. In the simplest case of trans- 
lation invariance an orthogonal Fourier transformation 
U = {Uki} diagonalizes all matrices, 



m k 1 ^ 



/./. ' 



m k ui k 5kk> 

Kn,k $kk' 



= E^ Uki Uk '3 A % ' 

= E« Uki Uk 'i B % > 



so that 



EE ln 

fc n=l 



v n + (! + m k K nyk )u;l 



(45) 
(46) 
(47) 

(48) 



the renormalization coefficients of Eqs. (36) and ( |37j ) be- 
come 

Ay = U kl U kj A fc , dj = U kl U kj a k , (49) 
where the variances of the k-th mode, 



X k 



a k 



m k 



E 



1 + m k K n , k 



E 



(3m k v\ + (1 + m k K n , k ) lo\ 



(50) 
(51) 



generalize those of Eqs. ( |l^ ) and (|2p|); for instance, the 
'on-site' renormalization coefficient D = Cu can be sim- 
ply expressed as D = (Cf) = N~ 1 J2k a k- The parti- 
tion function and thermal averages are to be evaluated by 
means of Eq. (|35|), where, of course, the effective poten- 
tial and the double-bracket average are to be understood 
as the LCA ones. Note that, since the LCA A's do not 
depend on q, averages of observables involving only mo- 
menta are trivially evaluated as (0(p)) = (0(p)\. 



IV. APPLICATIONS 
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A. Drude response model 

In the case of Ohmic dissipation the memory is Marko- 
vian, rj(t) = r]5(t — 0), corresponding to t)(z) = 77 = const 
and then to k„ oc |n|, so one can see that the coor- 
dinate fluctuation a(q), Eq. (|2fj|), is divergent, as well 
as /j(<t~), Eq. (18). This is due to the unphysical as- 
sumption of a vanishing response time from the dissi- 
pation bath: in the standard SPR model such diver- 
gence only affects the momentum fluctuation, 
leaving the possibility to meaningfully evaluate thermal 
averages of coordinate-dependent quantities ||] . As when 
considering momentum-dependent quantities R] in the 
standard case, in the anomalous SPR model we are there- 
fore forced to take into account the finite bath's response 
time. 

The Drude model of an exponentially decaying mem- 
ory, 



V(t) = 6(t)?7W D e" 



fj(z) 



?7 

w n + z 



(52) 



is the simplest one which is physically meaningful: the 
constant 77 characterizes the strength of the coupling with 
the dissipation bath, while the Drude frequency w D char- 
acterizes its bandwidth is the response time); the 
Ohmic limit occurs for lo^ — > 00. 



B. Single particle in the double-well 

Given a potential V{q), it is convenient to provide a di- 
mensionless formulation by identifying characteristic en- 
ergy and length scales, e and a. The first one could be 
a barrier height or a well depth, the latter is such that 
the change of V over its scale is comparable to e. The 
dimensionless coordinate is x = q/a and the dimension- 
less potential v(x) is such that V(ax) = ev(x). If x m is 
the absolute minimum of v(x) and v" = v"(x m ), the har- 
monic approximation (HA) of the system is characterized 
by the frequency u>q — y e v" /ma 2 ; the ratio between the 
HA quantum energy level splitting huio and the overall 
energy scale e defines the dimensionless coupling param- 
eter 



fJ 



hojo 

e 



h z v" 



meo~* 



(53) 



Therefore the system is weakly (strongly) 'quantum' 
when g is small (large) compared to 1. Defining the di- 
mensionless momentum p x = ap/h, such that = i, 
the nondissipative part of the Hamiltonian (0) reads 



1 ~2 

9 Px 



v(x) 



(54) 



In what follows energies are given in units of e, lengths 
in units of a, frequencies in units of loq, and so on; the 
dimensionless temperature is t = l/(e/3). 



As for the dissipation part, to make contact with the 
phenomcnological behavior, note that the asymptotic 
damping function from the classical limit of the pseudo- 
Langevin equation (|J) is V"fj(z) — muj 2 f)(z) and has the 
dimension of a frequency [it is indeed the counterpart 
of 7(2) /m of the standard SPR model], so for the Drude 
form ( |52[) we will deal with the corresponding dimension- 
less input parameters 



The dimensionless term mn n reads then 

Tin 

mn n = T] uj d 



wbstef D =0hu>j2=(g/2t)<j D . 
Using the dimensionless parameter 



(55) 



(56) 



fix) 



(3huj(x) g u(x) 
2 ~ 2t w 



(57) 



in the place of lv(x), from the definitions (|18|), (pL9|) , and 
201) we obtain 



n=l 

v"t 



^nf + (l + a n )p{x) 



\{x) 
a(x) 



9* 



E 

— -CK 
OO 

E 



(7m) 2 



^n) 2 + {l + a n )p{x) ' 
1 + a n 



2v"t ^ (7rn) 2 + (I + a„)/ 2 (x) 



(58) 
(59) 
(60) 



where A = (p 2 ,))- 

The example we choose here is the double- well poten- 
tial, 



v(x) = (1 - .t 2 ) 2 



(61) 



which has two degenerate symmetric minima in x m = ±1 
with v" = 8. From Eqs. (p% and (Ph we obtain 



v e g(x) = (1 — x ) — 3 a (x) + t n{x) , 
f{x)^^[ix 2 +Za{x)-l] . 



(62) 
(63) 



Eq. ( |63D has to be solved self-consistently with Eq. (J60|) : 
we did this numerically; exact reference data can be ob- 
tained only for the nondissipative limit 77 = 0, e.g., by 
numerically solving the Schrodinger equation. 

In Fig. [l] we report the shape of the coordinate prob- 
ability distribution V(x) — (5(x — x)), for selected val- 
ues of the dissipation strength r) and bandwidth tZ> D , at 
the coupling g = 2; this gives a ground state energy 
e = 0.69797, with the next excited level at e\ = 1.04637, 
above the barrier. When dissipation is switched on V(x) 
tends to go farther away from the classical distribution, 
V c ~ e~ v ^' t , due indeed to the enhanced fluctuations 
of the coordinate, and the two-peaked structure is even- 
tually lost. This happens by rising either 77 or uj d : note 
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the strong influence of the latter in increasing the fluctu- 
ations, which diverge in the Ohmic limit lu d — > oo. The 
figure also reports the opposite result for the standard 
dissipative model, obtained as in Ref. Q, which clearly 
shows the suppression of the coordinate fluctuations to- 
wards the classical distribution. 




FIG. 1. Configuration density V(x) — {S(x — x)) of the 
double-well quartic potential for g = 2, t = 1, and different 
values of the damping parameters fj = 1,2 (as indicated on the 
picture) and uj d — 1 (long-dashed lines) and 5 (short-dashed 
lines). The solid line is the nondissipative effective poten- 
tial result for fj = 0, that agrees with the corresponding exact 
result (filled circles); the dotted curve corresponds to the clas- 
sical limit, and the dash-dotted curve reports the result of the 
standard SPR model from Ref. ^| for an Ohmic dissipative 
strength 7 = y/fjniJo) = 5. 

Typical results found for the average potential energy 
v(t) = (v(x)) are displayed in Fig. In the nondissipa- 
tive case, comparison with the exact data shows that the 
effective potential gives very accurate results, in spite of 
the strong coupling. The curves for the anomalous dis- 
sipative system show an increase of the potential energy, 
corresponding to the higher fluctuations of the coordi- 
nate, while for the standard SPR model the opposite oc- 
curs. 



C. 4> chain 

The quantum </> 4 chain with standard dissipation was 
studied in Ref. (?]]; the effective potential approach was 
previously applied in the nondissipative case jl7],[l2| and 
checked through quantum Monte Carlo simulations . 
It can be viewed as the discretized version of a contin- 
uum nonlinear field theory, described by the (undamped) 
Lagrangian 



N 



c 

2iv 



tt 2 

2 (ft -ft-i) 2 + -s-v(qi) 



(64) 



where a is the chain spacing, c is the 'relativistic' veloc- 
ity, ft is the gap of the bare dispersion relation, G is a 
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FIG. 2. Average potential energy (v(x)) of the double-well 
quartic potential vs reduced temperature t, for g = 2 and 
different values of fj and lo d . Curves and symbols as in Fig. [l]. 



constant with the dimension (mass x length), and peri- 
odic boundary conditions ensure translation symmetry. 
The local nonlinear potential 



v(x) = (l-x 2 Y 



(65) 



has two wells in x m = ±1 [with v"(x m ) = 8], so that 
the classical <f) A chain has two degenerate translation- 
invariant minimum configurations, = 1} and {qi — 
— 1}, as well as relative minima connecting the two wells, 
the static 'kinks'. In the continuum limit ia — > z the 
kink configuration is indeed q(z) = ±tanh[f2(z — zq)/c], 
so that it is localized with a characteristic length c/Q, and 
its energy is e K = 2G£lc/3. 

The ratio between the characteristic frequency Q of 
the quasiharmonic excitations of the system and the en- 
ergy scale e K defines the quantum coupling parameter 
Q = Q/s K = 3/(2Gc); the discreteness of the chain is 
measured by the kink length in lattice units, R = c/(f2a) 
(R — > 00 in the continuum limit). Moreover, t = T/e K 
is the reduced temperature used from now on. The most 
interesting features appear when kinks are excited in the 
system; for instance, they cause a peak in the classical 
specific heat at t ~ 0.2. 

Using the dimensionless quantities just introduced, the 
Weyl symbol for the undamped </> 4 Hamiltonian can be 
written as 



H Q 2 R 



V(q) 



— Y 

2R ^ 



-ylft-ft-i) +— 



(66) 
(67) 



where the momenta are such that [qi,pj] = iSij, and 



Eq. (|5J) gives m^ 1 = 2Q 2 i?e K /3 = m" 1 . 

As for dissipation, N identical independent environ- 
mental baths coupled to each momentum are assumed, 
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giving a diagonal kernel matrix K n ,ij — <5y K n , so - 
gives K ntk = K, 
the Drude form 



Eq. © 

We take K n — \u n \ fj(\v n \) with 
for fj; the dimcnsionless quantity 
mkKn k — mn n = a n can eventually be written as in 
Eq. (|56|), so that dissipation is characterized by the di- 
mensionless parameters w,, = and fj — 3i]/(2QR). 

From Eqs. (ff^ ) and (|46] ) the LCA renormalized fre- 
quencies are found to be 



fl 2 [1 - 3D(t) + 4i? 2 sin 2 (fc/2)] 



(68) 



while the relevant renormalization coefficient D(t) = Cj, 
follows from Eqs. © and (||), 



0(f) = 



EE 



1 + a n 



Q 2 R 1 

3* ^V^iM 2 + (!+<*») /* 



(69) 



where / fc = /?w fe /2 = (Q/2t) [l - 3L»(t) + AR 2 sin 2 (fc/2)] ; 
these two self-consistent equations can be solved numeri- 
cally, and it appears that the approximation is meaning- 
ful when 3D <§; 1. Eventually, the LCA effective poten- 
tial can be written as the original one, just replacing the 
local interaction v(x) by 

VcS (x) = [1 - 30(f) - x 2 ] 2 + 6D 2 {t) + f/2(f) , (70) 
with 



16i? 1 \ \ 



k n—l 



(7rn) 2 + (l+a w )/ 2 (rr) 
(to) 2 



(71) 



For any parameter set (t, Q, R, fj, w D ) the self- 
consistent computation of D and hence of the last term 
of jl, which takes a negligible computer time using a con- 
tinuum termination of the n summation, completely de- 
termines the effective potential v e s{x). The problem is 
then ready for a numerical evaluation of Eq. ([55]), that 
gives the partition function (setting 0=1) and the ther- 
mal averages of observables. We employed the numerical 
transfer matrix technique |]l9| , that reduces the configu- 
ration integral for a one-dimensional array with nearest- 
neighbor interaction to a secular integral equation. The 
evaluation of the latter is implemented numerically using 
a discrete mesh for the values of each degree of freedom. 
Temperature scans over the region of interest were per- 
formed and several thermodynamic quantities were cal- 
culated, taking the value 'per site' for the extensive ones. 
In particular, the internal energy u(t) — f(t) — td t f(t) 
and the specific heat c(f) = d t u(t) = —td 2 f(t) follow 
from the free energy per site, f(t) = —N~ 1 t In Z(t) by 
numerical derivation, while Eq. ( J35[ ) was used to calcu- 
late the thermal average of the squared site coordinate 
(q 2 ) and of the local potential (v(qi)) , the square nearest- 
neighbor displacement ((ft — <?i-i) 2 ), and the square mo- 
mentum (p 2 )- The quantities reported in Figs. ||, ^, and f| 
are evaluated for fixed values of the kink length R = 5 
and of the quantum coupling Q = 0.2, which give fairly 
strong quantum effects. For comparison, we also report 
the classical result corresponding to Q = 0. Moreover, a 
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FIG. 3. Renormalization coefficient D{t) [Eq. (£>y|)J vs re- 
duced temperature t, for different values of the damping 
strength rj, at fixed coupling parameter Q — 0.2, kink length 
R — 5, and Drude cutoff frequency w D = 1. Solid line: 
fj — (undamped); short-dashed line: fj = 1; long-dashed 
line: r) — 2; dash-dotted line: fj — 5. Note that D(t) increases 
when fj increases. 



fixed representative value of the Drude cutoff frequency 
lu d = 1 is used, in order to analyze the dependence upon 
the dissipation strength fj. 

As we already remarked, the pure-quantum fluctua- 
tions of the coordinate are made stronger by the anoma- 
lous dissipative coupling. This appears in Fig. [|, where 
the pure-quantum renormalization coefficient D(t) — 
(£ 2 )) of Eq. ( |69| ) is reported for different values of fj, 
starting from the nondissipative value fj = 0. The tem- 
perature dependence of the coefficient D(t) affects the 
effective potential and the free energy, allowing us to de- 
scribe how the classical behavior of thermal averages is 
affected by both quanticity and dissipation. 

In Fig. |i| the temperature behavior of the mean-square 
fluctuations of the site coordinate (q 2 ) is reported. At 
t = 0, in the classical case the coordinate lies in the min- 
ima {q 2 = 1} of the potential (|67|), while in the quan- 
tum non-dissipative case the value at t — is smaller, 
(q 2 ) = 1 — 3D ~ 0.84: this corresponds to the minima 
of v e ft(q) and reflects the fact that quantum fluctuations 
make the configuration to climb the barrier rather than 
the steeper walls. Switching on the temperature enhances 
the same effect and (q 2 ) decreases at finite temperature, 
until when at t > 0.5 many kinks are excited in the sys- 
tem and the coordinate distribution begins to spread to- 
wards the walls, causing the subsequent increase from a 
minimum value ~ 0.56 at t ~ 0.45, eventually collaps- 
ing onto the classical curve. At variance with the case 
of standard dissipation |7|, when the damping strength 
f) is switched on, a further enhancement of the distri- 
bution spread occurs: this effect could be qualitatively 
interpreted as an 'effective' increase of the quantum cou- 
pling and allows (qf) to reach smaller values, at smaller 
temperature. 

On the other hand, the fluctuations of the momenta 
are quenched by dissipation, and the role of the damping 
effects is nonpredictable on a simple basis if one con- 
siders thermodynamic quantities where both coordinates 
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FIG. 4. Mean-square fluctuations of the coordinate opera- 
tor (g?) vs reduced temperature t. Parameters and lines as in 
Fig. H. The dotted line represents the classical result (Q=0). 



and momenta enter into play. One of these is the spe- 
cific heat c(t). Its nonlinear part, namely its total value 
minus the corresponding (dissipative) harmonic contribu- 
tion, 5c(t) = c(t) — Ch(t), is very sensitive to the nonlin- 
earity of the system, since its value is zero in a harmonic 
approximation. The fact that our approach retains all 
classical nonlinear features is crucial for getting Sc, whose 
temperature behavior is reported in Fig. B. Switching on 
and increasing the strength of the environmental cou- 
pling fj the curves can be seen to go farther apart from 
the classical one, which is a nontrivial result. A physical 
explanation is that anomalous dissipation increases the 
uncertainty in the kink positions, substantially making 
them longer and effectively raising their mutual interac- 
tion. 



V. CONCLUSIONS 

In this paper we have discussed the problem of anoma- 
lous dissipation, previously considered by Leggett Q and 
derived the effective potential formalism for an inter- 
acting many-body system with this kind of dissipation, 
which can arise from a microscopic description of the 
system-bath coupling, as done for instance in Ref. JhJ. 

Our approach allows one to reduce quantum mechan- 
ical thermodynamic calculations to a classical-like con- 
figuration integral, where both quantum and dissipative 
effects are included. In order to deal with many degrees 
of freedom the necessary low- coupling approximation has 
been introduced. The latter, if the system's symmetries 
are exploited, results in very simple expressions for the 
renormalization coefficients appearing in the theory. This 
is shown in detail for the case of translation symmetry. 

Applications of the framework have been made to the 
single particle in double-well potential and to a discrete 
(j) 4 one-dimensional field, whose strong nonlinearity yields 
kink excitations that play a dominant role. In both ex- 



0.3 - 
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t 

FIG. 5. Nonlinear contribution to the specific heat, defined 
as Sc(t) — c(t) — Cji(f) (times the kink length R), vs reduced 
temperature t. Parameters and lines as in Fig. H. The dotted 
line represents the classical result (Q=0). Increasing fj, 5c(t) 
goes farther away from the classical behavior, as if Q were 
raised 0. 

amples a Drude-like spectrum of the environmental cou- 
pling has been chosen for dissipation, and its influence 
on thermal quantities has been analyzed. 

The effective potential is a unique tool for dealing with 
such a system, since any approximate theory must re- 
tain the strong nonlinearity, which has a mainly classi- 
cal character, and this rules out conventional perturba- 
tive approaches. The method proved to be useful for 
the above model systems, so it is expected that it will 
find application in up-to-date physical problems. A pos- 
sible application of the formalism regards the effect of 
the blackbody electromagnetic field [|10|] onto the phase 
diagram of two-dimensional Josephson junction arrays. 
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APPENDIX A: STANDARD SPR MODEL 

Eliminating the bath variables from the equations of 
motion for the standard SPR model (Q) one gets || 
the Langevin equation (^|), where the coupling with 
the bath is described by the memory function j(t) — 
Q(t) m(,u)J cos(u>i t), whose Laplace transform is 

7(z) = J2 m ^i .a T ; ,2 > ( A1 ) 

and by the Gaussian random force F(t), whose time cor- 
relators can be expressed in terms of the memory function 
by spectral relations J| depending on the bath's temper- 
ature. 
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The quantum thermodynamics, after the bath has been 
'traced out' in the path- integral for ([[]), is described by 
the density matrix la] 



p(q",q') = V[q(u)} e~ s ^ u ^- s M u K (A2) 

Jq' 

~ m f(u) + V( q (u))} (A3) 



S[q(u)} = 



m 

T 



and S 1 given in Eq. (]3[); the Matsubara transform of the 
kernel k(u) reads 



(A4) 



and is thus given by k n = \v n \ j(ji/ n \). For a harmonic 
potential, V(q) = mui 2 q 2 /2, the evaluation of the path 
integral can be performed exactly and the density matrix 
turns out to be a Gaussian with partition function and 
variances as follows Hi: 



Z„ = 



<P 2 )h = 



_^ oo 



BTilj 11 -v 2 + uj 2 + k n /m 

°° 2,11 

m \ ~t uj + k n /m 
3 ^ vI+lo 2 + k n /m ' 

n— — oo " 

i °° i 

_ y - 

dm ^— ' v 2 + lo 2 + k„ lm 



(A5) 



In the nondissipative limit, k n — >0, the summations 
give the well-known results Z — 1/(2 sinh/), (q 2 ) — 
Ji- coth/ and (p 2 ) = SfS coth/, with / = dhu/2. 



APPENDIX B: ANOMALOUS SPR MODEL 



while the fluctuating force term, 



(B6) 



is a Gaussian random process, whose correlations arise 
from the assumption that the initial values {qe,Pi} corre- 
spond to the equilibrium at the temperature T = 3~ x for 
the isolated bath, namely (qtqt>) = 5u> 2r ^ ^ coth , 

<ptPt>) = 5u> r ^f L coth^p , =s2 f . _ 

The thermal density matrix for the harmonic oscillator, 
V(q) = muj 2 q 2 /2, can be evaluated reducing the problem 
to the standard-dissipation one, just by performing the 
canonical transformation q — ► — p/(mtu) and p —> muiq\ 
to get exactly the same form of Eq. (Q) it is sufficient 
to replace pg — > mujqg, qj> — > —pt/ (into) , an d — > 
(mui) 2 / (mguj 2 ). One finds, of course, Eqs. (A5) with the 
exchange of the expressions (q 2 ) H — > (mw) (p 2 ) H and 
(p 2 ) H — » (mui) 2 (q 2 ) H . Due to the above redefinition of 
the bath masses {m^}, the kernel (A4) takes the form 
£;„ = (mw) 2 k„, with 



1 



m e v 2 + uf 



so eventually one has 



(A 



Bhui v 2 + (1 + mK n )uJ 2 

oo o 

m \ - u> 

P ^ vl + (1 + mn n )u> 2 ' 



1_ ~ 



1 + mK r , 



mB ^— ' ^ + (l + mK„)a) 2 

n— — oo u v 



(B7) 

(B8) 
(B9) 
(BIO) 



From (|4|) we get the equations of motion 

m« + V'(g) = V-[m«w?&-7'(g)] l (Bl) 



m f (g £ +w 2 %) = V'(g) 



(B2) 



the second one has the general solution qz(t) — qj(t) + 
qf(t), with the retarded form of the particular solution, 

mecuj q P l {t) = V'{q{i)) 

dt' cos [w € (t-0] d v V'(q{t')) , (B3) 

o 

and the homogeneous solution 



q]>(t) = qe cos(u^i) 



m£UJi 



sin(w^) . (B4) 



Inserting this result in Eq. ( p3l| ) one gets Eq. (||), where 
V (t) = 0(t) V — cos(a;/t) , (B5) 



APPENDIX C: ANOMALOUS INFLUENCE 
ACTION 

For the bath part of Eq. ([IJ) the Feynman-Kac expres- 
sion in terms of the momentum can be used: the defini- 
tion of the influence action is therefore 



P -s,W - _ 



l-]lfv\pt(u)]e- s *[te*<-«»], (CI) 



where Z B — J^(2 sinhw^) 1 is the partition function of 
the isolated bath, and the bath action reads 



5 n = 



E 



i 



pj(u) +uj 2 (p e (u) -p(u))' 



(C2) 



mi 



Performing the Fourier transformation into Matsubara 
components p(u) = Yl n P n e WnU [and same for pi{u)} 
one gets 
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s * = E E i^»i 2 + ^» - p-\ 



(C3) 



(^+^ 2 )IWn| 2 + 



Z/2 + U}\ 



where the variable shift pt n = pg n — -jirhp p n has been 

performed. The path integral over pe(u) gives just the 
partition function Z B , so one is left with the anomalous 
influence action 



SAP(U)} = - E K n\Pn\ 2 , 



(C4) 



n— — oo 



i.e., Eq. (||), where k„ is again given by Eq. (B7), so that 
comparin g w ith the Laplace transform of the memory 
function (|B5|), 



n 



(*)=£ — 



mi z 2 + u>f 
one has K n = \v n \ fj(\v n \) , i.e. Eq. (@). 



(C5) 



APPENDIX D: EVALUATION OF p FOR ONE 
DEGREE OF FREEDOM 

Let us start from the definition (|l4| ) of the reduced 
density in terms of the system's trial action, given by 
Eq. @ with the trial potential (|l|), 



hS 



du 



%(quPu-PuC[u) 



Pu 



W 



mu 2 {q u -q) 



2, 



2m 2 

+^{qoP(3-poqp) + i[(qf3-qo)p- (pp-pa)q] , (Di) 



and of the dissipative action (g) . In the calculations q is 
to be regarded as a fixed parameter and we can omit to 
explicit the dependence of w and oj 2 on it. Performing 
the shift q u — > q+q u and using the Fourier representation 



2^ 



exp 



i mLU 2 z q^ 
(D2) 



one gets 



p {p,q-q)=(3muj 2 e-P w J J 



(D3) 



0H 

du 



-^{quPu-Puqu 



Pi . muj 2 ql . 2 

1 imuj zq u 

2m 2 



+^(ioP0-poq0) + i[(qp-qo) p - (pp-po)£\ ■ (D4) 

where £ = q — q. To eliminate the linear term we shift 
now g„ — > g„ + iz, 



2^ 



p {p,q;q) = (3muj 2 e 0W 
so that ^2 is the harmonic-oscillator action, 



(D5) 



ns 2 



du 



^{quPu-Puqu) + -~~ 



muj 2 ql^ 



+ 2(ioPfj-Poqp)+'i[(qf3-qo)p- (pp-po) (f-^)]. ( D6 ) 

Now, the point is t o get rid of the path integral ap- 
pearing in Eq. (D5), that represents the Weyl symbol 
for the density operator of the harmonic oscillator with 
the dissipative action S'Jpn]. However, the result for the 
harmonic oscillator with anomalous dissipation is known 
from Eqs. (jB§), ( pTpj ) and Q: 



27re^ er p2 / 2X T e~ q ' ' l 2a i 



(D7) 



where the partition function fro m (|B8|) is written as Z — 
/(phuj) and p given as in Eq. Q1S|). One gets therefore 



Po(p,q;q)=f3muj 2 e-e w 



dz 
2^ 



(D8) 



This last Gaussian convolution affects the coordinate 
part and results in a Gaussian for £ with the pure- 
quantum variance a — a T — l/(mu) 2 (3) reported in 
Eq. (pp|): what is subtracted is ju st th e classical contri- 
bution, i.e., the n=0 term in Eq. (BIO). The final result 
is just Eq. ©. 



APPENDIX E: EVALUATION OF p FOR MANY 
DEGREES OF FREEDOM 

We evaluate here the path integral ( |3l| ) with the trial 
action So [i.e., Eq.(|3^) with the 'potential' (^4])]; the fixed 
argument of w(q) and of the matrix B(q) is omitted in 
the following. First, the delta function is represented as 



^Iof^)=det((3B 2 )J J ^ F 



e iS§du *zB 2 {q u -q) 



(El) 



and the exponential can be incorporated into the action 
So; then, performing the shift q u — ► q u + q + iz and 
introducing the fluctuation variable £ = q — q one gets 



po(p, <?;<?) 



,-/3 W 



det C c J (2n 



e 2 c *Pha(P.€-**) . 



(E2) 



where C" 1 = (3B 2 and 



^p,q)=Jv[p u ,q u ] e -5[ P „,^iP,«l-s 1 K,^;9]-5 I b u ] ) 



(E3) 
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with 



ft = 



t p u A 2 p u + t q u B 2 q l 



(E4) 



The path integral (E3) corresponds to the density ma- 
trix of a harmonic system with anomalous dissipation, 
p HA (p,q). The canonical exchange q u — > p u and p u — > 
—q u has the effect of turning the dissipation into stan- 
dard, ^[pj — * S'JqJ; performing the same transforma- 
tion onto the external variables (p, q) and exchanging 
the matrices A 2 <-> B 2 , an exact correspondence with 
the standard dissipative harmonic system is established, 
p HA (p, q) — p HS (q, —p), and the known result for a har- 
monic system with standard dissipation |Q ) can be used, 
eventually getting 



Pha(P><?) = 



where 



det(/3AB) y/(2n) N det C* H ^ir)" det A H ' 

(E5) 



Mh = 



C„ = 



A„ = 



det(^ + V r 

,.2N 



71=1 " n 

-° ,,2 , 



n— — oo 



1 oo 

1 £ B 



1 



and 



B [A 2 + K n ] B 



(E6) 
(E7) 
(E8) 

(E9) 



Eventually, the Gaussian convolution in Eq. (E2) washes 
out the n—0 component of C, leaving 



po{p,q;q) 



1 

2^/3 



1 



det A 



yJ{lTx) N det C \J (2n) N det A ' 



(E10) 
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